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Nitrogen-vacancy ensembles promise an unprecedented combination of sensitivity and resolution 
for magnetic field sensing under ambient conditions. However, their full sensing potential is limited 
by the inhomogeneous degradation of their collective coherence, while dynamical-decoupling limits 
sensing to the detection of AC-fields. Here we propose using the recently demonstrated cavity- 
protection effect for room-temperature, ensemble-based DC metrology. We identify a parameter 
space where strong-coupling can persist at room temperature and where the dressed states are 
limited solely by the dephasing of the cavity mode and the ensemble longitudinal relaxation, showing 
aT Hz“ 2 DC-held sensitivities to be theoretically achievable with further augmentation through 
Fano-interference. 

PACS numbers: 76.30.Mi, 76.70.Hb, 42.50.Ct, 42.50.Pq 


Magnetic metrology schemes employing nitrogen va¬ 
cancy (NV) spin ensembles promise an exceptional com¬ 
bination of sensitivity and spatial resolution, due to the 
increased fluorescence yield and scaling of their 

read-out spin-projection noise [1], where N is the num¬ 
ber of spins in the ensemble. Due to the ease of their 
coherent control under ambient conditions and their bio¬ 
compatibility, NV centres have emerged as particularly 
attractive quantum-bits for exploring quantum informa¬ 
tion storage and transferral [2, 3], as well as sensors for 
studying nano-scale electromagnetic helds and temper¬ 
ature gradients in biological specimens [4-6]. However, 
there are severe limitations with increased ensemble size 
based on the cumulative degradation of ensemble coher¬ 
ence through both intrinsic and extrinsic sources of inho¬ 
mogeneity. Conventional techniques used to suppress de¬ 
coherences resulting from inhomogeneity limits the over¬ 
all detection bandwidth of NVs in such a way that makes 
DC-helds below ~ 1 /iT difficult to measure [7, 8]. 

Here we propose employing dressed spin-ensemble 
states for DC-magnetic field sensing. The studied scheme 
would involve monolithic microwave resonators with high 
room-temperature quality factors (Q’s) that would fa¬ 
cilitate optically detected magnetic resonance (ODMR). 
We explore the feasibility of obtaining strong spin¬ 
coupling at room temperature and the DC-sensitivity 
of the resulting dressed states, with further enhance¬ 
ment proposed through Fano-based interference with a 
second cavity mode. The sensing described here em¬ 
ploys spin-resonance fluorescence-dependent detection, 
however phase-contrast based detection (Ramsey-type 
schemes) are equally applicable. Room temperature 
ensemble-QED and macro-scale microwave cavities have 
been considerably studied [9-14], and the cavities de¬ 
scribed in these references are anticipated to provide a 
basis for the proposals outlined in this article. 

The dominant decoherence source for NVs stems 
from dipole-interactions with the surrounding Over- 
hauser field, in addition to the local variations of the 


crystal held resulting from free charges, strain and surface 
states [15-17]. As a result, the dispersion of NV Larmor 
frequencies for an ensemble is inhomogeneously broad¬ 
ened. Furthermore, ensembles are intrinsically more sen¬ 
sitive to the amplitude and frequency huctuations of con¬ 
trol helds which poses difficulties for wide-held metrology. 
This is further exasperated by the random distribution 
of four possible crystallographic NV orientations. Pro¬ 
tocols such as Hahn echo, CPMG-n or XY-n have been 
routinely used to extend coherence times, but result in 
optimised sensitivity only at bandwidths inversely pro¬ 
portional to half the systems total coherence time (U), 
and when phase-locked with an external AC-held to max¬ 
imise the measured phase contrast [18-20]. With rigor¬ 
ous noise cancellation, sensitivity in the pT Hz “2 range 
for AC-helds has been demonstrated using a Hahn-echo 
pulse scheme [21]. 

While there is much incentive in developing optimised 
quantum-control and noise-discrimination protocols, al¬ 
ternative approaches have been explored and demon¬ 
strated in the context of hole-burning [22], weakly- 
coupled cavities masing between iris = 0 —)■ — 1 states 
[23], strongly coupled ensemble systems [24, 25], and 
recently a proposal of combining both strong-coupling 
and hole-burning [26]. In particular, the suppression of 
inhomogeneous-related decoherence in strongly coupled 
systems, the cavity-protection effect, has been theoreti¬ 
cally outlined for NV spin ensembles [25] and has been 
recently demonstrated experimentally using a supercon¬ 
ducting strip-line resonator [27]. 

The intention here is to utilise the cavity protec¬ 
tion effect in the scheme described in the Fig.l(a,b). 
Cavity-protection ensures that the identical constituents 
of a spin ensemble resonant with the cavity transition 
(or dispersed within the cavity bandwidth) are isolated 
from the in-homogeneously dispersed subset, provided 
that the coupling strength exceeds the sum of the cav¬ 
ity bandwidth and the inhomogeneous dispersion width. 
In particular, the strong-coupling suppresses the in- 
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FIG. 1. (Colour online) (a) Schematic of the cavity-ensemble system, where the dressed states ’tp±a.re accessible for an incoherent 
probe field with rate v, and the resulting fluorescence is detectable through the same or similar aperture, (b) Expected anti¬ 
crossing of the ODMR lines when tuned into resonance with a microwave cavity mode, with the splitting being a function of the 
ensemble size N. Scattering of the dressed state can occur with nearby phase-mismatched cavity modes (see Fig.2(c,i)). Inset 
is the intended pulse-scheme, (c, d, i) Time-dependent decay of (ct^) and (a^a) for an ensemble size of 10^° with k = 3 kHz ( 
Q = 10®), 7 i = 250 Hz, 72 = 1 MHz, 1 Hz (corresponding to a mode volume of 2 x 10“® cm® at a resonant frequency 
(aJc = 3 GHz), (c, d, ii) Steady-state values (a^}ss and (a^d)ss plotted as a function of incoherent pump rate and temperature, 
starting from a thermally equilibrated state (the same parameters as in (i) are used). 


homogeneously dispersed tails of the spin-distribution, 
effectively turning them into dark-states [28] . Such a sys¬ 
tem provides the distinct advantage of suppressing detri¬ 
mental broadening caused by (coherent) control-field er¬ 
rors, retaining coherences similar to those of single NVs, 
while providing the increased fluorescence yield of an en¬ 
semble measurement when probed incoherently. 

NV spins are prepared in the ground state (rus = 0) 
using an incoherent optical pump with a rate v, which is 
applied for a duration that is long enough to ensure that 
the ground-state spin population is maximised and has 
reached a steady-state [29] . This is followed by a coherent 
microwave 7r-pulse with a strength tj and a frequency 
(jjp that is resonant with the cavity and a chosen spin 
transition (in Fig.1(b) the -1-1 state is chosen), which is 
tuned into resonance with the cavity mode by applying 
an external magnetic field. When resonant, the resulting 
dressed states ip± are formed from a collective spin-wave 
(Dicke state) which posses a splitting fi and a degree 
of inhomogeneous suppression that is dependent on the 
initial spin dispersion and number of in-phase spins. This 
is in turn primarily determined by the thermal bath and 
the initial coherent/incoherent pump rates. The dressed 
states are then read-out using a saturating incoherent 
pulse, destroying the collective spin-wave and resetting a 
majority of the spins into the ground-state [29]. 

In general, a collective spin-wave in an ensemble can 
be created and sustained at pump rates where the num¬ 
ber of induced excitations for a given duration does not 
exceed the number of ensemble constituents, assuming 
all spins are identically polarised prior to excitation. In 
the solid-state this is difficult at elevated temperatures 


due to the constant exchange of energy between the lat¬ 
tice and the spins, yet it is still possible depending on 
the ensemble constituent size, collective dephasing rate 
and, most importantly, the amplitude of the polarising 
excitation. This is difficult to theoretically model in the 
solid state, except at low temperatures where one can 
justifiably collectively treat the spin ensemble under the 
Holstein-Primakoff (HP) approximation. This approxi¬ 
mation allows for the formal definition of the spin oper¬ 
ators in terms of boson operators, under the assumption 
that the number of excitations within a collective spin 
wave is much less than the number of its constituents. 
To justifiably apply this approximation to estimate the 
system output at room temperature, we are limited to 
theoretically explore a parameter space where the inco¬ 
herent preparation pump rate does not exceed the cavity 
and ensemble decay rates, as well as account for the ther¬ 
mal distribution of the ground state spins following the 
incoherent preparation pump, before the polarising co¬ 
herent TT-pulse. For the sake of the following derivations, 
we also assume the the 7r-pulse is weak enough to ensure 
that the HP-approximation holds. 

The interaction of a cavity mode with an ensemble of 
]V spins can be described by a general Hamiltonian using 
Pauli spin operators, and the dipole and rotating wave 
approximations: 

^ N N 

n=u}ca^a+ - '^ojk^l+'^gkia'^a-a'^d^) ( 1 ) 

fc k 

where h = 1, uj^. and gk are the transition angular fre¬ 
quency and coupling strength of the spin. The dy¬ 
namics which account for the cavity decay rate k, spon- 
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taneous emission 71 (= l/Ti) and coupling to the ther¬ 
mal environment h, polarisation dephasing 72 (= I/T 2 ), 
and loss through incoherent pump rate v in the con¬ 
text of ODMR, can be numerically calculated using a 
master equation for a reduced density matrix: ps = 

^K(n)\.Ps\ ~\~ •^ 7 i(ft)[Ps] + •^72 [Ps] (see 

Appendix A). Following similar justihcations as in [30], 
this is expanded neglecting third and higher-order corre¬ 
lations to explore the temperature and incoherent pump 
rate parameter space where a HP-based description re¬ 
mains valid. The relevant parameters in this case are the 
steady-state values for the inversion operator and 

the intra-cavity photon number {a^a)ss prior to the ap¬ 
plication of the coherent 7 r-pulse. It is assumed that the 
HP condition is fulfilled as long as {a^a)ss is below the 
spin number N in the ground state prior to polarisation 
into the rus = - 1-1 state, so as to allow polarisation to oc¬ 
cur dominantly by the coherent microwave pump rather 
than through incoherent remnant thermal excitations. 

As shown in [31], the full set of coupled equations need 
not be solved, as the spin population will have reached 
a steady-state before the coherent microwave pump is 
reapplied. The simplihed derivation results in four cou¬ 
pled equations and a closed pair of steady-state expres¬ 
sions (see Appendix A). Fig.l(c,d) shows a summary of 
both the temperature and incoherent pump-dependent 
dynamics obtained for (d^) and (a^a) for an ensemble of 
10^° NVs, where the time-dependent dynamics are calcu¬ 
lated assuming a thermally distributed spin population 
at t < 0. The dynamics are presented to show the com¬ 
bined effect of temperature, dephasing and incoherent 
pump rate on the dressed state phase coherence, as well 
as highlighting a limit where v exceeds the ensemble de¬ 
cay rate and destroys the coherent dressed state. As long 
as u < { 71 , k}, the number of excitations accumulated at 
the systems steady-state prior to the application of the 
coherent 7 r-pulse will be low enough to justify a HP-based 
description at any temperature when preparing and pop¬ 
ulating the dressed states. Given the exceptional rate 
of 7 i for NV’s [32] such excitation rates correspond to 
powers in the pW-fW range depending on the system, 
while higher powers in the preparation stage will dimin¬ 
ish the initial ground-state spin population and phase 
coherence. The steady-state values of the inversion op¬ 
erator also highlights the dominant effect of temperature 
irrespective of v, where the fraction of excited spins is 
diminished to a thermally mixed state with increasing 
temperatures through the relation extracted from the fit 
of Fig.l(c,ii): (ct^(T))s, = tanh (j^)(d^(0))^s. 

Eq.l may now be described using normalised collec¬ 
tive operators with a temperature dependent collective 
coupling strength while applying the HP-approximation. 
Assuming all spins are equally coupled to a magnetic 
mode with their quantisation axis aligned with the cav¬ 
ity field, a cooperative coupling strength can be defined 
in terms of the temperature tanh-dependency extracted 


from the fit of Fig.l(c,ii) as a temperature-dependent 
(and therefore A^-dependent) coupling strength: flT,N = 
g\/N tanh wq/ 47 rfcf,T. Here g ex ^e \/is the 
individual magnetic coupling field strength defined in 
terms of the gyromagnetic ratio 7 e 2.8 MHz/G and 
the effective mode volume 14®. By accounting for the 
temperature-dependant number of spins, the disordered 
portion of the spin ensemble is excluded, thereby only 
accounting for the spins that occupy the same eigen¬ 
state. Eq.l is rewritten as: T-Lhp = oJc&'^d -I- + 

'l2k9k{b^d — ho)) [33], and using the input-output for¬ 
malism, the steady-state cavity output tc{io) can be de¬ 
fined when 7 i is smaller than the inhomogeneous width 
(71 < A en = 72 ) [28, 34] as: 

f «;/2 

^ Aw-I-i(K-I-7i-I-27ril|, j^p(a;, Ae„)) — ^A 

’ ( 2 ) 

where Aw(= — w) is the detuning between the in¬ 
put field and the cavity resonant frequency, p{u}, Ae„) 
is the spectral distribution with a full width at half 
maximum of Ae„, and A is the non-linear Lamb-shift. 
ODMR will depopulate all ground-state levels and high¬ 
light both strongly-coupled and uncoupled spin levels. 
The expected spectra is plotted in Fig.2(a) for a Gaus¬ 
sian spin distribution, highlighting the occurrence of both 
cavity-pulling and suppression of emission detuned from 
the dressed state energies when fully on-resonance (see 
Appendix B for expression) . These effects are respec¬ 
tively diminished and enhanced as the coupling strength 
is increased and the dressed-state manifestation occurs 
at larger detuning (a small coupling strength was used 
to highlight these effects in Fig.2(a)). The temperature- 
dependent cavity output as a function of ensemble size 
is plotted in Fig.2(b). The dressed-states linewidth can 
be extracted from Eq.2 at the dressed-state frequencies 
(w = flT,N) ss: 

_ K-l- 71 -|- 27 rHf,_^p(r 2 T.Ar, Ae„) 

— 2 

Eq.3 is plotted in Fig.2(c) for various spectral dis¬ 
tributions using the temperature-dependent coupling 
strength. Fp, and thus the coherence time of the dressed- 
state, is directly dependent on the spectral dispersion. 
In the context of single NVs where dephasing is homoge¬ 
neous (Lorentzian) the influence of coupling to a cavity is 
detrimental. This is in contrast to ensembles which suffer 
from inhomogeneous broadening (the larger/denser, the 
more Gaussian-like), for which strong coupling will sig¬ 
nificantly suppress 72 . In the limit of infinite coupling, 
Eq.3 is reduced to (k -I- 7 i)/ 2 , for which the onset can 
be seen in Fig.2(c). It is worth mentioning that even for 
mode volumes in the order of the resonant wavelength, 
the dipole-coupling between NV’s ( 7 g/ro/ 2 /i « 3.3x 10“^^ 
Hz cm^) will always be comparably negligible to the col¬ 
lective coupling strength. 
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FIG. 2. (Colour online) (a) Detuning-dependent {Au = u)c — uj) ODMR plot of a strongly coupled system with the axis divisions 
normalised in units of the inhomogeneous linewidth Ae„. (b) Surface plot of tc(tu) (Eq.2) as a function of ensemble size N 
for two different temperatures, highlighting the need for larger collective coupling strengths to generate dressed-state at higher 
temperature, (c) Eq.3 plotted for three different distributions (a 1:4 Psuedo-Voigt function is plotted), along with the resulting 
linewidth of a scattered dressed-state (Eq.4) assuming a Gaussian spin distribution. In this plot k = 3 kHz (Q = 10®), 71 = 250 
Hz, 72 = Aen = 1 MHz, g — 1 Hz, and T = 300 K. 


This limit can be further augmented through Fano- 
interference when considering scattering between the 
dressed state and another uncoupled resonant mode. 
Cobnsidering e.g. circular resonators, these possess a 
dispersion of whispering gallery modes whose phase will 
be mismatched depending on the relative mode num¬ 
ber difference and spatial distribution. Scattering of the 
coupled mode can occur with any other mode order, ei¬ 
ther in the same cavity, or with an additional externally 
coupled cavity. If a dressed state is tuned to the fre¬ 
quency of a neighbouring mode by varying the coupling 
strength, or evanescently interfered with a mode of an¬ 
other resonator, Fano-type scattering is expected to oc¬ 
cur, as well as magnetically induced transparency (MIT) 
when on-resonance [35]. An example of a detuning plot 
of a dressed state interfering with a 7r-mismatched mode 
is shown in Fig.3(a,i) with extracted ODMR spectrum 
plotted in Fig.3(a,ii). The resulting linewidth of the 
MIT is determined by the phase mismatch and the on- 
resonance overlap integral between the cavity mode and 
the dressed-state mode, and takes the approximate form: 


Fm/t 


FpF, 
Fp -b 


( 4 ) 


where (j) is the phase mismatch and Fg is the scatter¬ 
ing mode linewidth. This results in a spectral window 
which can possess a narrower linewidth then the dressed- 
state but which exhibits half its contrast, as plotted 
in Fig.3(a,ii). Alternatively, the asymmetric Fano line- 
shapes, which emerge with detuning, do not suffer from 
reduced contrast. The gradient of these slopes exceeds 
that of both the individual dressed state and the MIT 
window, as highlighted Fig.3(b). 

Both bare and Fano-scattered dressed-states can be 
applied in schemes where sensing is conveyed through 
measured fluorescence contrast (see Appendix D). Exter¬ 
nal magnetic fields can be inferred through the spin Lar- 
mor frequency shift in their presence via Avesr = leBo^ 
where Bq is the projection of the magnetic field along the 
principal NV quantisation axis. Given that the dressed- 
state frequency is explicitly dependent on the detuning 
between the cavity resonance and the central frequency of 
the spin distribution, it is fair to assume that the dressed- 
state will shift with an identical gyromangetic ratio to 
that of an uncoupled NV spin. The frequency dependent 
fluorescence can be quantihed in terms of the expected 
measured number of photons from an ensemble. During 
read-out, the ensemble is incoherently probed to satura- 
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FIG. 3. (Colour online) (a,i) Detnning-dependent (Acu = flT,N — o^scat. mode) ODMR plot of a dressed-state with a scattering 
mode, where the dressed-state’s cavity mode quality factor is equivalent to that of the scattering mode {Qp = Qs)- (a,ii) Plot of 
normalised extracted spectra from the coloured dotted lines in (a,i). (b) Plot of the maximum slope in the resulting spectra as 
a function of detuining, normalised to the the slope of the individual dressed-state. The maximum slope of the MIT linewidth 
is in the case worse then that of the un-scattered dressed state, while the parameters used are identical to those for the plots 
in Fig.2. 




FIG. 4. (Colour online) (a) Log-scale DC-sensitivities at 300 K as a function of Qp and a given density of NVs collectively 
coupling to a cavity mode with an effective mode volume of 1 cm^ (g « 0.04 Hz at Wc = 3 GHz). The dashed red line highlights 
a pessimistic room temperature limit for GHz resonator Q’s. The parameters used here are 71 = 250 Hz, 72 = Aens = 1 
MHz. (b) Log-scale DC-sensitivities at 300 K for the MIT probe defined in Eq.3 for the same parameters. The dashed red line 
highlights Qs = 10®. (c) Log-scale DC-sensitivities for Qp = 10®, measured using the point of maximum obtainable slope of an 
asymmetric Fano lineshape. 


tion, thereby generating emission from both coupled and 
uncoupled spins is collected. Approximating the point of 
maximum slope of the resulting spectra as a function of 
r, sensitivity can be defined through: 

1 r 

SBqdmr — — „ ^ (5) 

le CVC 

where C is the contrast of the peak height from the 
off-resonance intensity, and (( is the shot-noise limited 
detected photon number for which C = rjoetQe • ?7ext • 
N ■ tm ■ Icyc ■ Vnvqe, where detector quantum efficiency 
VDetQe (70 %), photon extraction efficiency (50 %), 
the number of photon emitters Af, the read-out pulse 
duration time tm (1 s), the total cyclic transition from 


excited triplet-states through the singlet-states to the 
triplet ground-state (m^ = 0) "fcyc 130 MHz), and 
the total quantum efficiency of the radiative decay rate 
from the m-s = 0 state tjnvqe (80 %) have been included. 

Fig.4 plots the predicted DC-held sensitivity for Fp 
(calculated using the coupled modes quality Qp), Tmit 
(calculated using Qp and a scattering mode quality Q^ 
= 10®), at 300 K, with an individual coupling strength 
calculated using a mode volume of 1 cm®. The predicted 
sensitivity spans from nano to below atto T Hz “2 but 
obtaining Qp > 10® at room temperature is technically 
demanding, beyond which stabilisation of the resonance 
becomes challenging. However, even considering a pes¬ 
simistic limit of Qp’s at 10^, sensitivities around aT Hz“ 5 
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seem to be achievable, given the demonstrated achievable 
densities of 10^® cm“® [36]. Using the MIT window, the 
requirements are lessened, as seen in Fig.4(b) for a Qs of 
10®. An almost linear independence of Qp is achieved up 
until Qp > Qs- Both plots highlight a decrease in sensi¬ 
tivity after a certain ensemble density before increasing 
further below aT which highlights the transition 

into the cavity-protected regime. Further log-scale im¬ 
provement is seen in Fig.4(c) when using the asymmetric 
Fano-lineshape when Qs > Qp = 10® for the same param¬ 
eters used in Fig.3(a), showing sensitivities approaching 
zT Hz-i 

In summary, the calculated dynamics highlight the 
manifestation of strong coupling and the diminishing 
influence of temperature for increasing NV ensemble 
sizes. Using a collective coupling strength and assum¬ 
ing that large ensembles will posses a dominantly Gaus¬ 
sian distribution, a temperature-dependent, steady-state, 
analytical expression for the cavity-protected dressed- 
state linewidth was derived, which was used to estimate 
ODMR based DC-held sensitivities. This predicted aT 
Hz “2 sensitivity for densities approaching 10 ^° cm“®, 
with further augmentation investigated through the pos¬ 
sibility of scattering with external resonant modes show¬ 
ing zT Hz “2 may be achievable. These results high¬ 
light the prospect of using both optical and spin degrees 
of freedom under ambient conditions to further explore 
cavity-QED and metrology schemes using NVs. 

We are grateful to Fedor Jelezko, Viatcheslav Do- 
brovitski, Dara McCutcheon and Jonas Neergaard- 
Nielsen for fruitful discussions and feedback. This work 
was supported by the Danish Innovation Foundation 
(EXMAD project). 

Appendix A: Expressions for and {a^a)ss 

The primary intention here is to justify the Heisenberg- 
Langevin solution for the steady-state system output us¬ 
ing the Holstein-Primakoff approximation, for a large 
{N > 10®) ensemble coupled to a single cavity mode at 
room temperature. This allows for an analytical expres¬ 
sion of the steady-state output , unattainable through 
direct numerical integration of the master equation or 
a Monte Carlo wave function approach, which accounts 
for system temperature and ensemble spectral dispersion 
shape (for N two-level systems, the density matrix will 
possess a dimensionality of [ 2 ^ x 2 ^] which becomes nu¬ 
merically intractable when N > 10® through both these 
procedures). 

The system is first explored using a reduced master 
equation to assess and obtain an expression for the ef¬ 
fect of temperature on the total inversion, as well as the 
steady-state number of remnant excitations, prior to pop¬ 
ulating the spin-ensemble states. It is assumed that the 
correlation time of the heat bath is much shorter then 


the system evolution time-scale as well as that the sys¬ 
tem and bath are initially uncorrelated. The Markov- 
approximation is therefore employed, in addition to as¬ 
suming that all ensemble spins are identical. The Liou- 
villian operators described in the main text are defined 
as: 

i^n{n)[Ps] =K{n -I- l){2dpsa^ - {a'l'a,/9s}) 

+Kn{2d^Psd — {dd\ Ps}) ( 6 ) 

1 ^ _ 

^ii(n)[Ps] =-[ii{n + l)^{2d^ Pscr+ - {ct+o-^ ,p 4 ) 

k 

N 

- {d-fc<T)^,pJ)] (7) 

k 

N 

^i2[Ps] =y Y^i^lPsdl - Ps) (8) 

k 

N 

^v[Ps] =--^(2(7+p,ct- - {&-&+, Ps}) (9) 

k 

Here the thermal excitation of the mode and ensemble 
is described through the temperature and frequency de¬ 
pendent average number of thermal photons: h(w,T) = 

huj 

l/(gfci,T _ 2 ). The operator expectation values are de¬ 
rived through cyclic permutation of the expectation val¬ 
ues and reduced density matrix trace Tr[Ops]. As the sys¬ 
tem is prepared and read-out using the incoherent probe, 
the total phase-invariance of the system can be assumed 
to remain unchanged for the instances considered here, 
namely the preparation stage before the coherent probe 
is applied. Under this assumption, and only considering 
correlations up to second order, the coupled rate equa¬ 
tion reduces down to a set of four coupled equations from 
which the inversion and cavity occupation number form 
a closed pair of equations: 

(d^) =-(71 -I- u)(( 7 ^)( 2 h -I- 1) -I- (71 - v) 

—i2g{{d'^d) — c.c.) ( 10 ) 

(aid) =2«:(h — (d’^d)) — igN{{a~a^) — c.c.) (11) 

(d-dt) =-i(wc + uik){d~a)) - ]^ig{{d^) + 1 ) 

-ig{N-l){d+d~)-ig{d^){a^d) (12) 

(d+d“) =—ig{(T^)[{a~ a)) — c.c.) 

-( 71 -bn)(d+d") (13) 

where the tilde annotation signified complex frequencies 
which incorporates all the loss rates through the complex 
component. The dynamics extracted from these coupled 
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equations directly highlight an ensembles collective po¬ 
larisation and coherence for a given set of initial tem¬ 
perature and incoherent pump conditions, through the 
manifestation, or lack of, Rabi oscillations. The steady- 
state expressions for the inversion and cavity occupation 
number can then be extracted from the coupled equations 
and expressed as: 


_ (u - 7 i) - (4K/A^)((al'a),„ - n) 
(u + 7 i)( 2 h+l) 


(14) 


{a^d)ss =n - 


^(71 + v) 

4k 

;< t "=)„(1 -f 2 h ) 


- 7i \ 
V + Jlj 


(15) 


where in denotes the initial inversion and cavity photon 
number prior to incoherent pumping. Considering an 
initial cavity system in thermal equilibrium, {a^)in = 0 
and {a)d)in = n. 


Appendix B: Total System Output 

Having outlined the effect of temperature and incoher¬ 
ent pump strength, the initial Hamiltonaina may be sim¬ 
plified within these approximations using expressions for 
the collective coupling strength, collective spin operators 
and the Holstein-Primakoff approximation: 

^ , /-TT . 

(7^ = 1 — Wh K 

k 
N 

= = 2b^b - 1 

k 

N 

^ ~ ^ ^ 9k 

k 

The resulting Hamiltonian is defined as: 

N 

'ftfiP — -h E ujipb + n{b^d — bd^) 

k 

The equations of motion for the expectation value of the 
two operators (a) and (6) can then be defined along with 
the input-output relations to derive an expression for the 
complex transmission of the cavity: 

(a) =*[H^p,(a)]-^(a) 

=- + t(wc - w)^ (a) -I- H(S) - ?7i„ vE2 (20) 

(S) =*[H^P,(S)]-|^(6) 

=-(^y + i(wfe - w)^ (6) - 6ffc(a) (21) 

ritr ={d)y/K/2 ( 22 ) 

(23) 


(16) 

(17) 

(18) 

(19) 


where rjtr , ? 7 re and rjin are the transmitted, refiected 
and injected external field amplitudes. Rearranging the 
steady-state expectation values for the total cavity trans¬ 
mission gives: 


tc(w) 


_ ?7»n(zK/2) _ 

0Jc-(^-J2k 9k/i^k-U}) 


(24) 


To account for the spectral dispersion, the number of cou¬ 
pled NV’s are described as a continuous spectral density 
distribution p{uj) with a full width at half maximum Ae„ 
that is equivalent to the collective pure dephasing rate 
which includes both homogeneous and inhomogeneous 
contributions. The spectral density is initially expressed 
as: 


N 

= (25) 

k 


This is redefined in Eq.24 and simplified through the fol¬ 
lowing substitution: 


1 

(Wfc - w) 



duj' 


S(uj' - LOk) 

(oj' — UJ — 171 ) 


(26) 


The integral in the denominator is then expanded using 
the Sokhotski-Plemlj theorem in the limit where 71 <C 





(w — w') 


(27) 





p{^' 

(w — Ijj') 




resulting in the transmission expression: 


(28) 


tc(w) 


IKp 

Aw -I- i{K -I- 7 i -I- 27rH2p(w, Ae„)) — 


(29) 


where Aw = cOc — co, p{<jj) is defined in terms of the inho¬ 
mogeneous broadening A^ns, and A is a principal value 
integral in Eq.28, interpreted as the resulting non-linear 
Lamb-shift, and the injected external field amplitude is 
set to unity. Einally assuming that all coupling strengths 
are equal, = g'^N, and may be further defined in 
terms of the temperature dependence of the number of 
spins occupying the same level. This way, the ordered 
subset of the steady-state thermally disordered spin en¬ 
semble can be accounted for, by defining the temperature 
dependence of {d^)ss- This dependence is derived ana¬ 
lytically from Eq.l4, or can be extracted through fitting: 

^T,N= 9\/WHfyus (30) 



where pq is the magnetic permeability, 7 e is the spin 
gyromagnetic ratio and VeS is the effective cavity mode 


Vtr —^re Vin 




volume. The resulting linewidth can be extracted from 
the transmission function by analysing the function poles 
which occur at w = flx.N and is expressed as: 

K + 7i + 2TT^'^pfp{n T,N 7 ^en) /Qo'v 

tp = X (di) 


The spectral dispersion may be defined using any distri¬ 
bution function; here we use Gaussian, Lorentzian and 
a pseudo-Voigt function to obtain analytical expressions 
that can be directly solved: 


PG 



(33) 


expressed via: 



where r is the radius and n is the refractive index of the 
disk resonator. 


Appendix D: Magnetic Field Sensing 

A sensors measurement sensitivity to any external 
field, in this case magnetic, is fundamentally defined 
through: 


Pl 


27r((Ae„/2)2 -1- Aw2) 


(34) 


5B 


min 



/Jetm max 


dX 

dvESR 


(40) 


Ppv = ttpL + (1 - '&)PG, 0 < I? < 1 (35) 


The total emission expected to be highlighted through 
ODMR will include both coupled and uncoupled spin 
states. This implies that the total transmission should 
be defined as the sum of the output of both steady state 
operators (a) and (6). The summed expression, plotted 
in Fig.2(a) in the main article is given as: 


Iodmr{^)= 


2 


tci^') F tensip^') 


(36) 


tc{<^) + 


(Aw + iK)Ay/tc{u)) 


“t.at 


(37) 


Appendix C: Fano-Scattering 


where a{t) is the standard deviation of the sensors signal 
as a function of time and dS/dB is the response of the sys¬ 
tem to a changing external magnetic field B. 5S{t) « ^/P 
is the corresponding shot noise of the total detected num¬ 
ber of photons per measurement f3 = Xt^. This is a gen¬ 
eral relation that can be applied to most measurement 
procedures with fairly intuitive and straightforward sub¬ 
stitutions. 

For real measurements, the detected light is dispersed 
with a particular line shape and intensity that is set by 
the coherence time, scattering cross-section, and the total 
quantum efficiency of both excitation and detection of 
the polarised spin states. The point of maximum slope 
of this dispersion is where the largest shift in fluorescence 
will occur for a given change in external field, and can 
be approximately defined through the full width at half 
maximum of a symmetric dispersion through 


Fano-interference can be summarised as the resultant 
spectrum of energy dissipating from an initial state to a 
final state through two competing paths that are energet¬ 
ically indistinguishable. This can be derived straightfor¬ 
wardly by defining the initial Hamiltonian for two differ¬ 
ent cavity modes with distinct coupling strengths to the 
same ensemble. In the scenario described here, this can 
be well-approximated by the overlap integral of compet¬ 
ing emission profiles (which in this case are the dressed- 
state spectrum and the uncoupled mode spectrum), while 
accounting for their relative phase. To first-order ap¬ 
proximation, the resulting spectrum can be expressed 
through: 


Ipanaiu}) = |tc(w) + tens(w) + 6 *'^Gcatt(w)P (38) 

where tscatt(w) would typically be defined through a 
Lorentzian profile, and 4> is the relative phase difference 
for two competing modes. Looking at e.g. a circular 
disk resonator and exclusively analysing 1®* order ra¬ 
dial modes, the relative phase mismatch of two differ¬ 
ent modes at frequencies c/Xdress and c/Xscatt can be 


max 


dX 

dvESR 


cx 

AvVjr 


(41) 


where Vp is a distribution function-dependent constant 
(for a Lorentzian and Gaussian lineshape this is equal 
to 4/3-\/3 and •\/e/81n2, respectively). This relation is 
used to assess the relative enhancement of the asymmet¬ 
ric Fano lineshapes in the main text. 
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